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Abstract 

In this paper, we analyse the Einstein and Einstein-Maxwell bil- 
liards for all spatially homogeneous cosmological models correspond- 
ing to 3 and 4 dimensional real unimodular Lie algebras and provide 
the list of those models which are chaotic in the Belinskii, Khalatnikov 
and Lifschitz (BKL) limit. Through the billiard picture, we confirm 
that, in D = 5 spacetime dimensions, chaos is present if off-diagonal 
metric elements are kept: the finite volume billiards can be identi- 
fied with the fundamental Weyl chambers of hyperbolic Kac-Moody 
algebras. The most generic cases bring in the same algebras as in the 
inhomogeneous case, but other algebras appear through special initial 
conditions. 
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1 Introduction 



It has recently been shown QjJ |2| that the classical dynamics of the spatial 
scale factors (and of the dilatons if any) of D-dimensional gravity coupled to 
p-forms can be described, in the vicinity of a spacelike singularity, as a billiard 
motion, i.e. as the relativistic motion of a ball inside a region of hyperbolic 
space bounded by hyperplanes on which it undergoes elastic bounces. These 
results generalize the work by Belinskii, Khalatnikov and Lifshitz (BKL) |3] 
that discovered the chaotic behaviour of the generic solution of the vacuum 
Einstein equation, in four dimensional spacetime, near such a singularity. 
These authors showed that, near a spacelike singularity, a decoupling of the 
spatial points occurs which hugely simplifies the dynamical equations of the 
spatial metric, in the sense that, asymptotically, they cease to be partial 
differential equations to simply become, at each spatial point, a set of second- 
order non linear ordinary differential equations with respect to time. 

These equations coincide with the dynamical equations of some spatially 
homogeneous cosmological models which possess some of the qualitative 
properties of more generic solutions. For D — 4, the spatially homoge- 
neous vacuum models that share the chaotic behaviour of the more general 
inhomogeneous solutions are labelled as Bianchi type IX and VIII; their ho- 
mogeneity groups are respectively SU(2) and SX(2,IR). In higher spacetime 
dimensions, i.e. for 5 < D < 10, one also knows chaotic spatially homoge- 
neous cosmological models but none of them is diagonal 0. In fact, diagonal 
models are too restrictive to be able to reproduce the general oscillatory be- 
haviour but, as shown e.g. in [S], chaos is restored when non-diagonal metric 
elements are taken into account. 

A systematic way to study the asymptotic behaviour of solutions of Ein- 
stein's equations in the neighborhood of a spacelike singularity and in any 
spacetime dimension is provided by the Hamiltonian approach developed in 
[S| in which gravity is coupled to a collection of p-forms. In this framework, 
a billiard description of the asymptotic evolution of the scale factors is nat- 
urally set up and chaos follows from the fact that the billiard's volume is 
finite. In many interesting examples, some of them related to supergravity 
models [HUE], the billiard can be identified with the fundamental Weyl cham- 
ber of a Kac-Moody algebra and the reflections against the billiard walls 
with the fundamental Weyl reflections which generate the Weyl group (or 
Coxeter group): accordingly, the finiteness of the volume, hence also chaos, 
rely upon the hyperbolic character of the underlying Kac-Moody algebra; 
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this last property is established through its Cartan matrix or equivalently in 
its Dynkin diagram. 

The purpose of this paper is to analyse the billiard evolution of spatially 
homogeneous non-diagonal cosmological models in D = 4 or 5 spacetime 
dimensions, in the Hamiltonian formalism mentionned hereabove. Since one 
knows that the full field content of the theory is important in the character- 
ization of the billiard, we compare the pure Einstein gravity construction to 
that of the coupled Einstein-Maxwell system. The gravitational models we 
are interested in are in a one-to-one correspondence with the real Lie alge- 
bras - a complete classification based on their structure constants exists for 
d = 3 and d — 4, (JD = d + 1)0 - and we restrict our analysis to the 
unimodular ones, because only for such models can the symmetries of the 
metric be prescribed at the level of the action. We proceed along the lines 
defined in but with a special concern about the role played by the con- 
straints. Indeed, while in the general inhomogeneous case, the constraints 
essentially assign limitations on the spatial gradients of the fields without 
having an influence on the generic form of the BKL Hamiltonian, in the 
present situation, they precisely relate the coefficients that control the walls 
in the potential. Consequently, the question arises whether they can enforce 
the disappearing of some (symmetry or electric or magnetic) walls as already 
do the vanishing structure constants with gravitational walls. They could 
prevent the generic oscillatory behaviour of the scale factors. The answer 
evidently depends on the Lie algebra considered and on its dimension: for 
example, while going from the Bianchi IX model in d = 3 to the correspond- 
ing U3S3 model in d — 4, the structure constants remain the same but the 
momentum constraints get less restrictive. Hence generic behaviour is easier 
to reach when more variables enter the relations. We find that, except for 
the Bianchi IX and VIII cases in D = 4, symmetry walls (hence off-diagonal 
elements) are needed to close the billiard table: thereby confirming, in the 
billiard picture, previous results about chaos restoration. Moreover, we find 
that when the billiard has a finite volume in hyperbolic space, it can again 
be identified with the fundamental Weyl chamber of one of the hyperbolic 
Kac-Moody algebras. In the most generic situation, these algebras coincide 
with those already relevant in the general inhomogeneous case. However, in 
special cases, new rank 3 or 4 simply laced algebras are exhibited. 

2 We will use the notations of MacCallum, we refer to ^2] for translation to other 
notations 
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The paper is organized as follows. We first adapt to the spatially homo- 
geneous case that part of the general Hamiltonian formalism set up in jU] 
necessary to understand how, at the BKL limit, the billiard walls arise in the 
potential. We explicitly write down the form of the momentum constraints in 
the generalized Iwasawa variables and analyse their meaning for each of the 
3 and 4 real unimodular Lie algebras as well as their impact on the billiard's 
shape. For the finite volume billiards, we compute the scalar products of 
the gradients of the dominant walls using the metric defined by the kinetic 
energy and show that the matrix 



is the generalized Cartan matrix of an hyperbolic Kac-Moody algebra. 

2 General setting 

2.1 Spatially homogeneous models, Hamiltonian 

In this paper, we are specially interested in d = 3 and d = 4 dimensional 
spatially homogeneous models equipped with a homogeneity group simply 
transitively acting; these models are known to be in a one-to-one correspon- 
dence with the 3 and 4 dimensional real Lie algebras and have been com- 
pletely classified |T21 1T3] . We restrict our analysis to the unimodular algebras, 
i.e. those whose adjoint representation is traceless 3 , that is C l ik = , since 
only for these homogeneous models do the equations of motion follow from a 
reduced Hamiltonian action in which the symmetry of the metric is enforced 
before taking variationnal derivatives. 

We work in a pseudo-Gaussian gauge defined by vanishing shift N l = 
and assume the D = d + 1 dimensional spacetime metric of the form 



where x° is the time coordinate, t is the proper time, dt = —Ndx°, and N 
is the lapse. For definiteness and in agreement with the choice made in j^J, 
we will assume that the spatial singularity occurs in the past, for t — 0. 

3 The group Adjoint representation is unimodular. 





ds 
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The gravitational dynamical variables are the components of the d 
dimensional spatial metric in the time-independent co-frame {ou l = u l jdx^ 
invariant under the group transformations 

du^—C'^Au"; (2.2) 

the C^'s are the group structure constants. The metric gij(x°) depends only 
on time and may contain off-diagonal elements. With use of g = detgij, one 
defines the rescaled lapse as N = N/ ^fg. 

In the spatially homogeneous Einstein-Maxwell system, there is besides 
the metric, an electromagnetic 1-form potential A and its 2-form field strength 
F = dA. In the temporal gauge A = 0, the potential reduces to 

A = Aj uj j . (2.3) 

In the Hamiltonian framework, we assume the potential itself to be spatially 
homogeneous 4 so that its space components in the oji frame are functions of 
x° only: Aj = Aj(x°). Accordingly, its field strength takes the special form 

F = dA = d Aj dx° Ad/-- Ai C l jk cu j A u k (2.4) 

which shows the links between the components of the magnetic field and the 
structure constants. Hence, from the Jacobi identity, one infers that 

F 1{J C\ £] = 0. (2.5) 

The first order action for the homogeneous Einstein-Maxwell system can 
be obtained from the D dimensional Hilbert-Einstein action in ADM form 
after space integration has been carried out; this operation brings in a con- 
stant space volume factor that will be ignored hereafter. The action is given 
by 

S[ 9ij , ir ij , Aj, 7r j ] = J cfe°(7r% + •• '•I - NH). (2.6) 

4 This is more restrictive than requiring spatial homogeneity of the field strength; in 
the present analysis the difference only arises with regard to the magnetic walls which are 
always subdominant. 
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The Hamiltonian NH reads as 



H = K + M (2.7) 
K = tt% - ^jTrVi + (2.8) 

M = (2.9) 

where i? is the spatial curvature scalar defined, in the unimodular cases, by 
the following combination of structure constants and metric coefficients 

R = ~\ (C llk Cjik + ~ C^ k C ijk ), (2.10) 

where 

C m = g u C l jk and C^ k = g^ g km C\ m . (2.11) 

The equations of motion are obtained by varying the action (|2.fij) with re- 
spect to the spatial metric components gij, the spatial 1-form components Aj 
and their respective conjugate momenta tt^ and tc j . The dynamical variables 
still obey the following constraints: 

H pa (Hamiltonian constraint) (2-12) 
Hi = —C- , ik -K k j + TC : 'FijpaO (momentum constraints); (2.13) 

notice that the Gauss law for the electric field is identically satisfied on 
account of the unimodularity condition. 



2.2 Generalized Iwasawa variables 

In order to develop the billiard analysis, it is necessary to change the vari- 
ables, i.e. to replace the metric components g^ by the new variables f3 a and 
jV", defined through the Iwasawa matrix decomposition 

g = M T A 2 M (2.14) 

where M is an upper triangular matrix with l's on the diagonal and A is a 
diagonal matrix with positive entries parametrized as 

A=exp(-/3), f3 = diag(f3\f3 2 } ...,f3 d ). (2.15) 
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The explicit form of (J2.14|) reads 



Hii^e-WMW*. (2.16) 

0=1 

The /3's are often referee! to as the scale factors although they more precisely 
describe their logarithms. The 7V"'s measure the strenght of the off-diagonal 
metric components and define how to pass from the invariant {u/} co-frame 
to the Iwasawa co-frame {O a } in which the metric is purely diagonal 

6 a =N° ( J. (2.17) 

In this basis, one has for the components of the 1-form 

Aj = AcM a j- (2.18) 

The changes of variables (|2.16|) and f)2. 18j) are continued to the momenta as 
canonical point transformations in the standard way via 

tt% + -'A, = -„ r + Y, vi *M a j + £ a A a - (2.19) 

a<j 

In this expression, V{ denotes the momentum conjugated to TV" and is de- 
fined for a < j, S a denotes the momentum conjugated to A a . The Iwasawa 
components of the electric and magnetic fields are given by 

£ a = N a ^ , T ah = F l3 M\M\ (2.20) 

where M\ denotes the element on line-j, column-a of the inverse matrix 
M~ x . This matrix enters the definition the vectorial frame {e a } dual to the 
co-frame {0 a } by 

e a = X ] M{. (2.21) 

While shifting to the Iwasawa basis and co-basis, the structure constants of 
the group, which also define the Lie brackets of the vectorial frame {X{\ dual 
to the invariant co-frame {uj 1 } 

[X h X 3 } = -X k C%, (2.22) 

transform as the components of a Q)-tensor so that 

ie h e c ] = -e a C\, with d\ c = N\N\N k c C^ k . (2.23) 
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2.3 Splitting of the Hamiltonian 

Following jH], we next split the Hamiltonian H into two parts: the first one, 
denoted by Ho, is the kinetic term for the local scale factors f3 a ; the second 
one, denoted by V, is the potential and contains all the other contributions. 
Thus 

H = H + V (2.24) 

and 

H = ^G ab ir a ir b . (2.25) 
The total potential naturally splits into 



V = V s + V G + VA + V™ agn (2.26) 



where 

a<b 

is quadratic in the Vs and as such related to the kinetic energy of the off- 
diagonal metric components; one refers to it as to the "symmetry" potential. 
Vg vanishes in the case of pure diagonal gij. Next comes the gravitational or 
curvature potential 

V G = -9R=\ e~ 2 ^ d C't C'\ c + \ e -*r+*+V {C '% c f ) (2.28) 

a,b,c 

involving the structure constants in the Iwasawa basis defined by (J2.23|) . 
The last two terms in the potential correspond to the electric and magnetic 
energy: 

V el = ^£ a £ a = ^e- 2ea (£ a ) 2 (2.29) 



ymagn = 1 e -2£f =1 /3 C JT ft jrab = 1 e -2m ab ^ 

where, with the notations of [6J, 



e a = P a , m ab = ( 2 - 3 °) 

c^{ab} 
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2.4 BKL Limit, billiard walls 



As the above formulae explicitly show, i) the total potential exhibits the 
general form 

V = V((3,Ar,V 7 £,F) = J2 ca{M,V,E,F) exp(-2w A {{3)) (2.31) 

A 

where wa{/3) = WAb/3 b are linear forms in the scale factors and ii) apart 
from the first term in the right-hand side of the curvature potential (which 
will introduce subdominant walls), the prefactors ca are all given by the 
square of a real polynomial, implying ca > 0. As explained in jH], with the 
following gauge choice, N = p 2 = —f3 a [3 a = —p 2r y a 'J a , and in the BKL limit 
corresponding to p — > oo, the exponentials terms in the potential NV become 
sharp walls and may be replaced by 

lim [cap 2 e- 2 ^«] = e oo (-2 Wj4 (7)) (2.32) 

p—>oo 

where 9 ro is defined through 

e ~M={°cor>°o < 2 - 33 > 

and has the property of being invariant under multiplication by a positive 
factor. Accordingly, in this limit, the positive c^'s can be absorbed by the 
Gqc and the null ones can simply be dropped out; once this has been done, 
the Hamiltonian no longer depends on the variables x 

X e{Af,V,£,F} (2.34) 

that enter these coefficients. This also means that the BKL Hamiltonian 
only contains the scale factors (3 a and their conjugate momenta ir a hence, as 
can be immediately infered from their equations of motion, all the x's must 
asymptotically tend to constants. These constants are arbitrary untill and 
unless the constraints enter the play. 

Another way to see that, before the constraints are taken into account, 
the asymptotic values of the x' s are arbitrary constants is to consider the 
map 

: Xo - X°° = <Kxo) (2.35) 
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which sends the initial values Xo on their asymptotic BKL values, for a given 
set of initial data {/3o,7r }, and to use the fact that, at least locally, this 
map must be invertible. This statement is clearly true for an Hamiltonian 
which does not at all depend on the x' s - It is reasonable to assume that 
this property remains valid in the more difficult situation that we face here, 
where the x' s only asymptotically get out of the Hamiltonian. 

In the general inhomogeneous case, the constraints are essentially condi- 
tions on the space gradients of the dynamical variables and, as such, they 
introduce no limitation on the generic BKL Hamiltonian: that's why gener- 
ically the Hamiltonian contains all the walls. The same happens of course 
when the constraints are absent. In the spatial homogeneity context, how- 
ever, they have to be taken into account since they can influence the billiard's 
shape. Of course, the prefactors controling the presence of the curvature 
walls crucially depend on the homogeneity group under consideration: the 
more the group "looks" abelian, the less curvature walls are present. The 
momentum constraints further establish linear relations between the other 
wall coefficients c^'s in which all together the structure constants, the x' s 
and even the variables {j3, 7r} are mixed up; so the question arises whether, 
asymptotically, they are equivalent to the condition that some of the ca's 
vanish forcing the corresponding walls to disappear. That is what we will 
systematically investigate in the following. 

3 d = 3 homogeneous models 
3.1 Iwasawa variables 

In spatial dimension d = 3, using the simplified notations 



N\ = m, N\ = ri2, A/" 2 3 = n 3 (3.1) 
V\ = p u V\=p 2 , V\=p 3 , (3.2) 



the prefactors (V^Afj) 2 of the possible symmetry walls e 2 ^ ^ a - ) , b > a, in 
dZ2B , read 



for a = 1, b = 2 
for a = 1, b = 3 
for a = 2, b = 3 




(3.3) 
(3.4) 
(3.5) 



C23 = {Pzf- 
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The Iwasawa decomposition (|2.14|) provides explicitly 

9n = e- 2 ?\ gn = n ie ~^\ g 13 = n 2 e~ 2 ^ (3.6) 
922 = n\e-^+e-^\ g 23 = n^e'^ + n 3 e^ 2 (3.7) 

and its canonical extension reads 

27T 11 = -(7Ti + 2-riiPi + 2n 2 p 2 )e 2(31 - (n 2 7r 2 + 2n\n 3 p 3 - 2n 1 n 2 p 3 )e w2 

— {n\ + n\n\ — 2n 1 n 2 n 3 )7r 3 e 2 ^ 3 (3.9) 

2tt 12 = Pl e 2pl + (nivr 2 + 2nin 3 p 3 - n 2 p 3 )e 2p2 

+ (nml - n 2 n 3 )7r 3 e 2 P 3 (3.10) 

2tt 13 = p^-n^ + h-nma)^ (3.11) 

2tt 22 = -(7i 2 + 2n 3 p 3 )e 2 ^ 2 -n 2 3 7r 3 e 2p3 (3.12) 

2tt 23 = p 3 e 2 ? 2 + n 3 n 3 e 2 ? 3 (3.13) 

2tt 33 = -vr 3 e 2/33 . (3.14) 

In order to easily translate the constraints in terms of the Iwasawa variables, 
we also mention the following usefull formulae 

2tt 2 = pi (3.15) 

2tt 3 = p 2 (3.16) 

2tt 3 = n lP2 +p 3 (3.17) 

2n\ = n 1 (n 2 -n 1 ) + (e- 2 ^ 2 '^-n 2 1 )p 1 + (n 3 e- 2 ^-^-n 1 n 2 ^ 

+ (nin 3 -ra 2 )p 3 (3.18) 

27Tg = n 2 (7r 3 - 7ri) +n 1 n 3 {n 2 - tt 3 ) + [n 3 e~ 2(p - n x n 2 \p x 

+ [e -2(/3 3 -^) + rfte-WP-fi 1 ) - n l] p2 

+ [n^l-n^-n^- 2 ^-? 2 ^ (3.19) 

2tt 2 3 = -n 3 (n 2 -n 3 )+n 2Pl + (e- 2 ^-^ -nl)p 3 (3.20) 

3.2 d = 3 pure gravity billiards 

These spatially homogeneous models are known in the literature as the class- 
A Bianchi models; they are classified according to their real, unimodular, 
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isometry Lie algebra [3 EI] • 

1. Bianchi-type I: Cf- = Q,Vi,j,k. 

This is the abelian algebra. There is no spatial curvature and the con- 
straints are identically verified. Accordingly, the billiard walls are only made 
of symmetry walls f3 l — ,i > j , among which the two dominant ones are 
w 32 — (3 3 — (3 2 and W21 = f3 2 — (3 1 . This is the infinite volume non-diagonal 
Kasner billiard. Its projection on the Poincare disc is represented by the 
shaded area in figure 1. 

2. Bianchi-type II: C\ 3 = 1. 

This case is particularly simple because the constraints are easy to anal- 
yse. Indeed, the momentum constraints read 

7r\ = and ir\ = p 1 = and p 2 = 0. (3.21) 

That means, referring to ()3.3|) . ()3.4j) and ()3.5|) . that cu = and C13 = and 
that they clearly eliminate the symmetry walls W21 = (3 2 — f3 l and = 
/3 3 - p 1 - the last P3 remains free so the symmetry wall W32 = /3 3 — ft 2 is 
present. Moreover, the only non zero structure constants being C'l 3 = 1, one 
single curvature wall survives which is 2/9 1 . Since two walls can never close 
the billiard, its volume is also infinite. Its projection on the Poincare disc is 
given by the shaded area in figure 1. 




Bianchi I Bianchi II 



Figure 1: Bianchi I and II billiards are the shaded areas limited by curved 
lines or curvature walls and right lines or symmetry walls. 

3. Bianchi-type VI : C\ 3 = 1 = C 2 13 . 
The momentum constraints read 

7r 3 2 = , tc\ = and tt\ + tt\ = 0. (3.22) 
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They are equivalent to 



= p 2 (3.23) 
= p 3 (3.24) 
= ni ( V r 2 - vri ) + ( e - 2 (/ 32 -/ 31 ) - n \ + 2) Pl . (3.25) 

The first two, according to (|3.4j) and (|3.5|) . clearly tell that the symmetry 
walls = j3 3 — j3 2 and w^i = j3 3 — (3 l are absent from the potential. 
With ()3.24|) put in ()3.3|) . one sees that the coefficient of the third symmetry 
wall W21 becomes c\2 = G°i) 2 and the question arises whether the constraint 
f!3.25|) implies pf 3 = 0? In order to answer this question, we will exhibit 
an asymptotic solution for which is a constant 7^ 0. Hence, can be 
different from zero and the symmetry wall W21 = (3 2 — (3 1 that it multiplies 
is generically present. 

Let us build such a solution. We know that after a finite number of colli- 
sions against the curvature walls, the ball will never meet them again. At 
such a time, one can find a solution of Kasner's type with ordered exponents. 
Let the diagonal Kasner solution be 5 

3 

ds 2 = -dt 2 + t 2qi {dx 1 ) 2 (3.26) 
i=i 

this is an asymptotic solution of the Bianchi VI type if qi > and q2 > 0. On 
this diagonal asymptotic solution, we next perform a linear transformation 
which generates another asymptotic solution which is no longer diagonal: let 

G = L T G Kasner L] (3.27) 

G is the new metric, Gxasner is the diagonal metric given above. We select 
L in the form of a Lorentz boost 6 in the plane (1, 2) such as to produce but 
a single off-diagonal element, namely gi2'- 



5 Here, the Kasner exponents will be denoted qi in order to avoid confusion with pi, the 
momenta conjugate to n% 

6 L is not really restricted here to define a Lie algebra automorphism, that is, the 
structure constants need not be conserved because, at the limit considered, the curva- 
ture becomes negligible and remains negligible after the transformation; however the L- 
transformation here chosen conserves the Bianchi VI structure constants and is a gauge 
transformation |l(Jj. 
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cosh a sinh a 

sinh a cosh a | . (3.28) 


Accordingly, the new metric elements are 

gn = t 2qi cosh 2 a + t 2q2 sinh 2 a (3.29) 
gi2 = (t 2qi +t 292 ) sinh a cosh a (3.30) 

the other stay unchanged. From these expressions and ()3.6|) . ()3.7|) . ()3.8|) . one 
immediately extracts the t-dependence of n\ and of the /3's which explicitly 
gives 

#12 (t 291 + t 292 ) cosh a sinh a . „ „ 1 . 

^i = — = o ^ 9 — (3.31) 

Sli t 2 *? 1 cosh'' a + t 2 ^ sinh 2 a 

and 

2G9 2_^i) = (t 2qi cosh 2 q + t 2q2 sinh 2 a) 2 

In order to discover the time behaviour of pi, we make use of the equation 
of motion for n\ which, since p\ just appears in the Hamiltonian through the 
coefficient Ci 2 , simply reads 

n 1 = ^ = ^^= Pl e- 2 ^ 1 \ (3.33) 

GST GEE GET 

remembering that [Bj dt — —y/gdr and that the sum of the Kasner ex- 
ponents is 1 such that y/g = t. Hence, the momentum pi is given by 
pi = —rii e 2 ^ 2- ^ 1 ); assuming q% < one finds that, in the vicinity of t = 0, 
it behaves as 

pi ~ C + B t 2{q2 ~ qi) + ... (3.34) 

where C = pf is not constrained to be zero. This result allows us to state 
that, in general, the symmetry wall 1021 = f3 2 — (3 1 is present, beside the two 
curvature walls 2/3 1 and 2/3 2 ; the dominant walls being w 2 ± = (3 2 — (3 1 and 
2/3 1 , they don't close the billiard. 

Notice that, for t — > 0, the exponential in the constraint ()3.25|) goes to 
zero like 

e -2(/3 2 -/3 1 )^ cosh -4 at 2( 92 - (?1 ) ) (3.35) 

as announced by the BKL analyzis when the coefficient of the wall c\i 7^ 
and also that the limit t — > of the product p\ exp — 2{(3 2 — (3 1 ) exists and 
is zero. Consequently, (J3.25)) admits a clear limit when t — *■ 0. 



13 



m 

Bianchi VI (generic) Bianchi VI (part.) 

Figure 2: Bianchi VI Billiards 

Had one taken advantage of the gauge freedom to assign the initial values 
n\ and p\ such as nf = 0, then the constraint (j3.25j) would have imposed 
the last symmetry wall to be absent and the billiard, here limited by the two 
curvature walls, to be open. 

These examples show that the presence/absence of symmetry walls may 
depend on initial conditions or on gauge conditions; nevertheless, in the cases 
mentionned hereabove, a billiard with an infinite volume remains a billiard 
with an infinite volume. 

4. Bianchi-type VII: C\ 3 = 1, C 2 13 = -1. 

This case is very similar to the preceding one. The only changes are i) 
that the third constraint in (J3.22)) has to be replaced by (-7T 2 — 7i\ = 0) and 
ii) that its translation into the Iwasawa variables in (J3.23)) now reads 

= m(7r 2 - tti) + (e- 2 ^ 1 ) - n\ - 2) Pl . (3.36) 

Accordingly, we can apply the same reasoning as before and conclude that 
generically, ^ and that the symmetry wall w%i = j3 2 — (3 1 appears in 
the potential. 

5. Bianchi-type IX: C\ 3 = 1, C 2 31 = 1, C\ 2 = 1. 

This case and the next one deserve a particular treatment because the 
structure constants are such that all curvature walls, namely 2/3 1 , 2/3 2 , 2/3 3 , 
appear and these three gravitational walls already form a finite 7 volume 
billiard. 

7 This situation is very specific to the homogeneous models in D = 4; in higher space- 
time dimensions, the number of curvature walls allowed by the structure constants is not 
sufficient to produce a finite volume billiard. 
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Let us first discuss the generic case. The momentum constraints take the 
form 

7T 3 2 - 7T 2 3 = 0, ir\ - 7r 3 = 0, vr 2 - tt\ = 0; (3.37) 
and in terms of the Iwasawa variables, they become 

n 3 (n 2 - 7r 3 ) - n 2 pi + riip 2 

+ ^ e -^-P 2 ) + n l + 1]p3 = (3.38) 



n 2 {-n 3 - tti) + nin 3 (7r 2 - 7r 3 ) + [n 3 e 2{/3 ' pl) - n 1 n 2 ]pi 

+ [n^ 2 - n 2 n 3 - nie- 2{/33 ~ /32) ]p 3 = (3.39) 



nifa - ttO + [e" 2 ^ 1 ) - n 2 - 2] Pl 
+ [ra 3 e~ 2(/3 ) - nin 2 ]p 2 

+ [nin 3 - n 2 b 3 = 0. (3.40) 

Remember that we already know that the billiard has a finite volume, hence 
the question is no longer to state between chaos or non chaos but rather to 
define more precisely the shape of the billiard. We cannot copy the reasoning 
made for Bianchi VI since we have here to account for an infinite number of 
collisions. 

In order to study the implications of the constraints (|3.38|) - (|3.4(Jj) . we 
shall rely on the heuristic estimates made in jH], where the asymptotic be- 
haviour of the variables is analysed in the BKL limit. From that analysis, it 
follows that, when p — > oo: i) the 7r a 's go to zero as powers of 1/p, ii) the n^s 
and the p^s tend to constants n°° and p°° up to additive terms which also go 
to zero as powers of 1/p and iii) that the exponentials either vanish (if the 
walls are present) or oscillate between zero and oo (if they are absent). What 
we infer from this information is that if the pf were not strictly zero then 
the constraint system hereabove would have no BKL limit. This is obviously 
wrong because a constraint need be obeyed all the time; so the limit must 
exist. Now we know that the p^s go to zero, can we say something more 
about the walls that they multiply? Compared to the exponential growth of 
the wall, the power decrease of its prefactor is not fast enough to prevent the 
symmetry wall to appear in the Hamiltonian. 
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Accordingly, the billiard's edge is formed by the leading symmetry walls 
u>32 and W21 and by the curvature wall 2/3 1 . On the Poincare disc, it is 
one of the six small triangles included in the larger one bordered by the 
curvature walls. Its Cartan matrix is that of the hyperbolic Kac-Moody 
algebra AE 3 = A^ A 

2-10 
-1 2 -2 
0-2 2 

already relevant in the general inhomogeneous case. Its Dynkin diagram is 
given in figure 3. 

o — 00 

Figure 3: Dynkin diagrams of the A± A Kac-Moody algebra 

Other interesting cases exist with less symmetry walls, which require spe- 
cific initial conditions. 




B. IX (generic) B. IX diagonal B. IX number 2 




Figure 4: Bianchi IX Billiards 



1. No symmetry wall at all. This situation is the one mentionned above; 
it happens when the metric is assumed diagonal always, hence for n\ = 
0,n 2 = 0,n 3 = 0; the solution of the momentum constraints is then 
Pi = 0,p 2 = 0,p 3 = 0. This assumption is consistent with the equations 
of motion. The Cartan matrix of this billiard is given by 

2 -2 -2 \ 

-2 2 -2 (3.42) 
-2-2 2 / 

and its corresponding Dynkin diagram is number 1 in figure 5; the 
associated Kac-Moody algebra is hyperbolic, it has number 7 in the 
enumeration provided in reference [T4] . 
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2. One symmetry wall. This happens when one chooses the initial data 
such that rii = 0, n 2 = and pi = 0, p 2 = 0; then, according to the 
equations of motion, these variables remains zero all the time. The 
billiard is closed by two curvature walls 2/? 1 and 2f3 2 and the symmetry 
wall W32- On the Poincare disc, its volume is half of that of the triangle 
made of curvature walls. The Cartan matrix is given by 

2 -2 \ 

2 -2 (3.43) 
-2-2 2 J 

and its Dynkin diagram is number 2 in figure 5. It characterizes the 
third rank 3 Lorentzian Kac-Moody algebra in the classification given 
in P3J. 




1 2 



Figure 5: Dynkin diagrams of special algebras met in Bianchi IX models 

These last two Kac-Moody algebras are subalgebras of A AA [TT] . 

Remark: let us recall that the reflections on the faces of the billiard gener- 
ate a Coxeter group which is identified with the Weyl group of the associated 
Kac-Moody algebra. The larger the set of walls, the larger the number of gen- 
erators. Since in the non generic cases the set of walls is a subset of the one in 
the generic case, the associate Coxeter group is a subgroup of the generic one. 

6. Bianchi-type VIII: C\ 3 = 1, C 2 31 = 1, C\ 2 = -1. 

The analysis of this case follows closely the previous one. The sign change 
in the structure constants modifies the constraints as follows 

7T 3 2 + 7r 2 3 = 0, 7^+7^ = 0, 7^-7^ = (3.44) 

and induces some sign changes in their Iwasawa counterparts. The conclu- 
sions are those of the Bianchi IX model: the generic billiard is the one of the 
algebra A AA . 
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3.3 d = 3 Einstein-Maxwell billiards 

The momentum constraints (j2.13|) of the Einstein-Maxwell homogeneous 
models generally (except for the abelian Bianchi I) mix gravitational and 
one-form variables; in comparison with the pure gravity case, i) no constraint 
remains which clearly forces the prefactor of a symmetry wall to be zero and 
ii) additional terms of the type itmF^ 1 ' enter in the constraints system. 

Accordingly, generically, besides the curvature walls allowed by the struc- 
ture constants, one expects all symmetry, electric and magnetic walls to be 
present. The dominant ones are the symmetry walls W21 = (3 2 — (3 l ,ws2 = 
j3 3 — (3 2 and the electric wall e\ = (3 1 which replaces the curvature wall of the 
pure gravity case. They close the billiard whose Cartan matrix is 



The Dynkin diagram is dispayed in figure 6; the corresponding Kac-Moody 
algebra is the hyperbolic A^ A algebra which is the Lorentzian extension of 
the twisted affine algebra also encountered in jjj. 



4 d = 4 homogeneous models 

Our billiard analysis of the four dimensional spatially homogeneous cosmo- 
logical models has been carried out along the same lines as for d = 3: results 
relative to pure gravity models and those relative to the Einstein-Maxwell 
homogeneous system will be given separately. 




(3.45) 




18 



4.1 Extension of the notations and Iwasawa variables 



In spatial dimension d — 4, we introduce the matrix variables 



( 1 


"12 


"13 


n u 


\ 





1 


"23 


n 2 A 










1 


"34 












1 


/ 



(4.1) 



The Iwasawa decomposition of the metric extends beyond the 3-dimensional 
formulae listed in (|3.6p . ()3.7|) . (j3.8j) . through the additional components 

-2/3 1 



#14 
#24 
#34 
#44 



riu e 

n 12 n u e~ 2/31 + n 24 e" 



2/3 2 



7113 "-14 e 



-2/3 1 



+ "23 "24 e 



14' 



-2/3 2 



+ 71 



34 1 



+ "34 e 
+ e 



2/3 3 
-2/3 4 



(4.2) 
(4.3) 
(4.4) 
(4.5) 



The momentum conjugate to (3 a is written before; the momentum 

conjugate to Afj = n a j is denoted V{ = p> a and is only defined for j > a. 

The prefactors c a b = (V{Afj) 2 of the symmetry walls e~ 2 ^ b > a, in the 
Hamiltonian are explicitly given by 



21 +P 31 "23+P 41 "24) 2 



Cl2 = (P 

C13 = (P 31 +P 41 "34) 2 

C14 = (P 41 ) 2 

C23 = (P 32 +P 42 "34) 2 

C24 = (P 42 ) 2 

C34 = (P 43 ) 2 . 



(4.6) 

(4.7) 
(4.8) 
(4.9) 
(4.10) 
(4.11) 



Once the change of dynamical variables has been continued in a canonical 
point transformation, the momentum contraints still express linear relations 
among the tt^s which translate into linear relations on the momenta ir a ,p ia ; 
their coefficients are polynomials in the n i's times exponentials of the type 
e -2(/3 6 -/3 a )^ with b > a. The explicit form of the constraints depends on the 
model considered. 



4.2 d = 4 Pure gravity 

We label the various 4 dimensional spatially homogeneous models according 
to the classification of the 4 dimensional real unimodular Lie algebras given 
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by M. MacCallum [12]: they are 



1. Class Z71[l, 1, 1]: C\ 4 = A, C\ = fi, C\ A = i/,with A + // + i/ = 0. 
One can still set A = 1 except in the abelian case where A = fi = v = 0. 

2. Class U1[Z, Z, 1]: C\ 4 = — |, C 2 14 = l,C\ 4 = — 1, C 2 24 = — |, 
C\ A — \i. n = is a special case. 

3. Class [A [2, 1]: C\ 4 = C 2 14 = 1, C 2 24 = -§, C 3 34 = fi and /i is 
or 1. 



4. 


Class C/l[3]: 


r 2 - 1 r 3 

O 14 — 1) <^ 24 — 


: 1. 






5. 


Class Z73J0: 


b 23 ~ ^ b 31 " 


' J-j ^ 12 - 


: -1. 




6. 


Class U3I2: 


^4 i 

b 23 ~ — X ' ° 31 


- 1 r 3 

— L > ° 12 


= 1. 




7. 


Class 17351 


or rf(2)®«(l): 


^23 = 1) 


C 2 3l = 1, 


C l2 = -1 


8. 


Class U3S3 


or sw(2) © m(1): 


C*23 = 1) 


c\ x = 1, 


C 3 2 = 1. 



Let us mention that for all of the four dimensional algebras, except of 
course the abelian one, the structure constant C'\ A 7^ and consequently 
that the curvature wall 2/3 1 + (3 2 is always present. 

Our analysis leads to the conclusion that, from the billiard point of view, 
the previous models can be collected into two sets: the first one has an open 
billiard, the second one has a finite volume billiard whose Cartan matrix 
is that of the hyperbolic Kac-Moody algebra A exactly as in the general 
inhomogeneous situation. 

The first set contains the abelian algebra and Ul[l, 1, lj^o.-ij Ul[2, 1] 
and U1[Z, Z, lj^oi ah the others belong to the second set. Because explicit 
developments soon become lenghty and since the reasonings always rest on 
similar arguments, we have chosen not to review systematically all cases as 
for d = 3 but rather to illustrate the results on examples taken in both sets: 

1. As a representative of the first set, we take U1[Z, Z, lj^o- 
The momentum constraints read 

tt\ = vr 4 2 = fJ,1T\ = -|(7ri + 7T 2 2 - 2 7T 3 3 ) + 1t\ 7T 2 = 0; (4.12) 
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in terms of the Iwasawa variables, they become 



p 41 =p 42 = //p 43 







(4.13) 



and 



- - (27T 3 - 7Tl - 7T 2 ) + ni 2 (7T 2 - 7Ti) + 

+ (e- 2 ^ 2 ^ 1 ) - n 2 2 - 1) p 21 + (n 23 e-W 2 ~^ - n 12 n 13 + ^ 



3n 23 /i 32 
+ (ni 2 n 23 - rii 3 H — ) p 



0. 



y 1 

(4.14) 



The first three constraints ()4.13|) clearly indicate that the symmetry walls 
(3 4 — (3 3 , j3 A — /3 2 , /3 4 — (3 1 are absent from the potential. What information 
can one extract from the fourth constraint (|4.14j) ? In order to answer this 
question, we will proceed as for the Bianch VI model: we will exhibit an 



asymptotic solution with prefactors c 23 = {p%) and ci 2 = (p 2 ^ + n^p 3 ^) 
different from zero. We here again face a situation in which the billiard 
motion is made of a finite number of collisions against the curvature walls. 
After the last bounce, the motion can be described by a Kasner type solution 



ds 2 



-dt 2 + t29i ( dxi f 



(4.15) 



i=i 



which is a suitable asymptotic solution at the conditions qi > 0, q 2 > and 
qs > (no condition on q^ is required because there is no curvature wall 
involving /5 4 ). We next perform a linear transformation 



G = L T G 



Kasner-' 



(4.16) 



in order to get another asymptotic solution which is non-diagonal. Here, for 
simplicity, L can be taken in the form that produces the off-diagonal elements 
we are intersted in, 

/ h h 

m\ m 2 
r\ r 2 




h 







1 1 



(4.17) 



Then using the equations of motion for ni 2 , rii 3 and n 23 we obtain 



P 21 + n 23 p 31 



n 12 e 



2(/3 2 -/3 1 ) 



(4.18) 
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and 



V = n 2 3 e 



2(/3 3 -/? 2 ) 



(4.19) 



By direct application of the formulae given in reference [H], section (4.2), 
assuming the following order qi < q 2 < qs for the Kasner exponents, we 
find that for t — > 0, (j4.18J) and ()4. 19|) tend to constants that depend on the 
parametrization of L but that will generically be different from zero. 

Consequently, our conclusion is that, in the generic case, the dominant 
walls of the billiard are the symmetry walls W32 = P 3 — P 2 , W21 = /3 2 — (3 1 
and the curvature wall 2/3 1 + /3 2 ; it is indeed an open billiard. 

Notice that in all cases of the first set, the non vanishing structure con- 
stants are all of the form C'\ h with a,b = 1, 2, 3. It is easy to check that 
this forbids all curvature wall containing f3 A and therefore that the remaining 
curvature walls cannot be expected to close the billiard. 

2. As a first representative of the second set, we take U1[Z, Z, l]^=o 

The constraints are given by ()4.12|) . (j4.13J) and (j4.14|) for // = and 
one immediately sees that, compared to the preceding case, one constraint 
drops out leaving the symmetry wall w i3 = /3 4 — j3 3 in place. To conclude 
the analysis of the last constraint, we can also provide a solution with non 
vanishing prefactors for the symmetry walls. Since the structure constants 
do not play an important role in the explicit construction hereabove, the 
solution obtained with /i ^ can also be used for fi = if one assigns 
particular initial condition such that p 43 = and ri 34 = 0. Again, we end up 
here with the following list of dominant walls: W43 = /3 4 — /3 3 , W32 = /3 3 — /3 2 , 
W21 = (3 2 — (3 l and 2/3 1 + (3 2 . The Cartan matrix is that of the algebra AQ A 
as previously announced. 



3. Another interesting example of the second set is provided by U3S3 

Its homogeneity group is the direct product SU(2) x U(l) so that this 
model appears as the four dimensional trivial extension of the Bianchi IX 
model: since the structure constants are the same as in the 3-dimensional 




o 



Figure 7: Dynkin diagram of the A algebra 
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case, the momentum constraints do not change, when expressed in terms of 
the metric components and their momenta 

7ri-7r? = 0, 4-7r? = 0, 7rf-7r? = 0; (4.20) 

their number does not change either but, in terms of Iwasawa variables, they 
involve much more terms than their d = 3 counterpart given in ()3.38|) , (|3.39|) 
and flUD} 

- "13(^1 - Ks) + n 12 n 23 (n 2 - tt 3 ) - (n 12 n 13 - n 23 e~ 2{ -^~ p ^p 21 

- (n? 3 - r^e" 2 ^-^) - e"- 2 ^ 1 ) - l)p 31 
+ (^12^3 - n 13 n 23 - n 12 e~ 2{l3i ~ p2) )p i2 

- (n 13 n 14 - n 23 e- 2 ^-^)n 24 - naie'^V 

+ (^12^23^24 - "-14^23 - ni 2 n 34 e ^ P P ')p 

- (n 12 n 23 n u - n l2 n 2i - n 13 n u + n M )p 43 = (4.21) 



- (ni 2 7ii 3 - n 2 3e~ 2 ^ 2 ~ pl) )p 31 + {n 12 n 23 - n 13 )p 32 

- (n 12 n u - n 24 e- 2 ^ 2 - 0l) )p A1 + (n 12 n 24 - n 14 )p 42 = (4.22) 



- n 23 (n 2 - vr 3 ) + n 13 p* - n 12 p 31 - (n 2 3 - e' 2 ^'^ + l)p 32 

(n 23 n 24 - n u e- 2 ^ 3 -^)p A2 + (n 23 n u - n 2i )p A3 = 0. (4.23) 

These constraints are of course i) linear in the momenta, ii) the only exponen- 
tials which enter these expressions are build of w 32 = /3 3 — P 2 , w 3 i = /3 3 — (3 l 
and w 2 i = /3 2 —/3 1 and hi) as before, their coefficients are exactly given by the 
square root of the corresponding c^'s in the potential, namely y/cu, y/ci3, and 
y/o^. Accordingly, the question arises whether these equations are equiva- 
lent to C12 = 0, c 2 3 = and c\ 3 = 0. If, asymptotically, the exponentials go 
to zero and can be dropped out of the constraints, the remaining equations 
can be solved for p 2 , p ,p 32 in terms of the other variables among which 
figure now, not only the 7Tj — tcj already present in the 3-dimensional case 
which asymptotically go to zero, but also the asymptotic values of the p Al1 s 
which remain unconstrained. It follows that none of the solutions of the 
above system is generically forced to vanish so that all the symmetry walls 
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are expected to be present. The absence of a wall can only happen in non 
generic situations with very peculiar initial conditions. 

We expect this result to become the rule in higher dimensions for trivial 
extensions like SU(2) x SU(2) x ... x £7(1) x ... x £7(1); the billiard will then 
become that of the Kac-Moody algebra A^ A relevant in the general Einstein 
theory. 

We can nevertheless provide a particular solution with initial data obeying 
ni2 = ni3 = nu = and p 21 = p 31 = p 41 = 0. These values are conserved 
in the time evolution. In this case, the leading walls are the symmetry walls 
W43, W32 and the curvature ones 2j3 2 + /3 3 , 2/? 1 + (3 2 . The Cartan matrix is 



/ 2 



V 



-1 
2 

-1 
-1 





-1 
-1 
2 



\ 



/ 



(4.24) 



The billiard is characterized by the rank 4 Lorentzian Kac-Moody algebra 
which bears number 2 in the classification given in ^3]; its Dynkin diagram 
is drawn in figure 8. 



Figure 8: Dynkin diagram of algebra number 2 in the subset of rank 4 
The D = 5 results are summarized in the following table: 



£71[1,1,1]^ =A=0 

£71[l,l,l]^ 0> -i 
£71 [Z, Z, 1]^ 


non chaotic 


£71[1,1,1] M= - 1 


chaotic 


£71 [Z, Z, 1] M=0 




£71 [2,1] 




£71 [3] 




£73/0 




£73/2 




£7353 




£7351 
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Chaos or non chaos for D = 5 models 



4.3 d = 4 Einstein-Maxwell models 

As in the three dimensional case, because of the presence of the electromag- 
netic field in the expression of the momentum constraint (J2.13j) . no symmetry 
wall can be eliminated. Moreover, the electric walls always prevail over the 
curvature ones if any. The billiard is in all cases characterized by the fol- 
lowing set of dominant walls = (3 A — /3 3 , w^2 = P 3 — P 2 , W21 = f3 2 — (3 1 
and e\ = (3 1 . Its Cartan matrix is that of the hyperbolic Kac- Moody algebra 

o — o — 0^0 

Figure 9: Dynkin diagram of the GQ A algebra 



In the general analysis of the billiards attached to coupled gravity + p- 
forms systems, one could assume 2p < d without loss of generality, because 
the complete set of walls is invariant under electric-magnetic duality. This 
invariance may however not remain in some spatially homogeneous cases due 
to the vanishing of some structure constants which lead to incomplete sets of 
walls so that other, a priori unexpected Kac-Moody algebras, might appear. 
An illustrative and interesting example is given in D = 5 by the Einstein + 
2-form system governed by the abelian algebra. Here, the dominant walls 
are i) the symmetry walls /3 A — (3 3 , (3 3 — f3 2 and j3 2 — /3 1 and ii) the electric 
wall (3 1 + j3 2 . This billiard brings in the new hyperbolic Kac-Moody algebra 
carrying number 20 in [14] . whose Dynkin diagram is given in figure 10 
hereafter. 



o 



Figure 10: Dynkin diagram of the algebra relevant for the Einstein + 2-form 

system 
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5 Conclusions 

In this paper, we have analyzed the Einstein and Einstein-Maxwell billiards 
for all the spatially homogeneous cosmological models in 3 and 4 dimensions. 
In the billiard picture, we confirm that in spacetime dimensions 5 < D < 10, 
diagonal models are not rich enough to produce the never ending oscillatory 
behaviour of the generic solution of Einstein's equations and that chaos is 
restored when off-diagonal metric elements are kept. Chaotic models are 
characterized by a finite volume billiard which can be identified with the 
fundamental Weyl chamber of an hyperbolic Kac-Moody algebra: in the 
most generic chaotic situation, the algebra coincides with the one already 
relevant in the inhomogeneous case; this remains true after the addition 
of an generic homogeneous electromagnetic field. Other algebras can also 
appear for special initial data or gauge choices: in fact, these are all the 
simply-laced known hyperbolic Kac-Moody algebras of ranks 3 and 4, except 
the one which has number 3 among those of rank 4 listed in [T3j. The billiard 
of a non chaotic model is even not a simplex. 
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